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Supplementary  Materials 

Text  SI:  Transcriptional  Signatures  of  Inflammation  and  Immune  System  Detected  in  Early 
Period 

Text  S2:  Transcriptional  Signatures  Associated  With  Injured  Muscle  Tissue  Extracellular  Matrix 

Text  S3:  Injured  Muscle  Tissue  Microenviromnent  Regulates  Existing  and  Invading  Cellular 
Migration,  Proliferation,  Differentiation  and  Phenotype 

Text  S4:  Activation  of  Muscle  Repair  Machinery 

Text  S5:  Late  Transcriptional  Programs  Activated  in  Response  to  Traumatic  LLMI 
Fig.  SI.  Immuno-histological  analysis  of  muscle  after  cryo-injury. 

Fig.  S2.  a)  Example  of  reproducibility  from  RNA-Seq  datasets.  Biological  replicates  of  in  vivo 
isolated  RNA  samples  from  24  hours  after  injury.  Sequencing  libraries  were  prepared  from  the 
isolated  RNA  and  a  strong  correlation  (R2  =  0.95)  was  observed  for  two  representative  examples, 
b)  Validation  of  RNA-Seq  data  with  qPCR  comparing  the  injured  versus  uninjured  samples  for 
the  3  hour  time  point,  c)  Validation  of  RNA-Seq  data  with  qPCR  comparing  the  injured  versus 
uninjured  samples  for  the  10  hour  time  point,  d)  Validation  of  RNA-Seq  data  with  qPCR 
comparing  the  injured  versus  uninjured  samples  for  the  24  hour  time  point. 

Fig.  S3.  Number  of  differentially  expressed  genes. 

Fig.  S4.  Gene  Sets  Activated  Early  After  Injury  (3-24  hrs)  a)  Enriched  KEGG  pathways  from 
differentially  expressed  genes  for  the  early  time  points  (3-24  hrs).  The  size  of  the  circle 
corresponds  to  the  number  of  significant  genes  with  each  enriched  pathway,  b)  Heatmap  of  genes 
associated  metabolism  expressed  in  the  early  time  period,  c)  Gene  expression  profile  of  Arginase 
1  (Argl),  which  is  a  molecule  emitted  by  type  2  macrophages  to  balance  inflammatory  reactive- 
oxygen  species  present  in  the  injured  tissue.  Red-Injured,  Blue-Uninjured 

Fig.  S5.  Examples  of  alternative  splicing  in  immune  network  of  genes. 

Fig.  S6.  Genes  associated  with  hematopoietic  activation,  proliferation  and  anti-inflammatory 
cytokines  activated  in  early  &  middle  periods. 

Fig.  S7.  Gene  sets  associated  with  extra-cellular  matrix  remodeling  &  satellite  cell  activation. 


Fig.  S8.  a)  Enriched  KEGG  pathways  from  differentially  expressed  genes  for  the  late  time 
points.  The  size  of  the  circle  corresponds  to  the  number  of  significant  genes  with  each  enriched 
pathway,  b)  Heatmaps  of  gene  sets  activated  in  late  period  categorized  by  function. 

Fig.  S9.  Area  plots  of  FDR  values  for  Gene  Annotation  (GO)  clusters  over  time.  GO  terms  were 
separated  into  major  functional  categories,  and  the  IQR  was  computed  for  the  tenns  in  each 
group.  The  grey  lines  indicate  the  25th  and  75th  percentiles,  red  lines  denote  the  25%  of  gene 
FDR  values,  the  blue  line  denotes  the  median  FDR  value,  and  the  green  line  denotes  the  75%  for 
FDR. 

Fig.  S 10.  Pathway-level  scoring  results  for  test  datasets  from  PCA  and  the  timepoint  signatures 
methods,  a).  Principal  component  analysis,  b)  Normalized  time  point  signatures  method  applied 
to  mean  pathway  expression  values.  Timepoint  signatures  for  6  test  samples  and  a  truth  training 
sample.  The  green  curve  indicates  the  control  sample  that  was  misclassified  by  the  algorithm  at 
the  pathway  level.  Time  point  signatures  for  the  9  time  points  (blue  lines)  as  well  as  6  test 
samples  (red  and  green  lines). 

Fig.  S 1 1 .  Differentially  expressed  gene  pathways  in  injured  samples  over  time. 

Table  SI.  Significant  GO  clusters  (FDR  <0.05)  identified  for  differentially  expressed  genes 
through  the  DAVID  annotation  tool.  Each  tab  represents  the  annotations  derived  for  each 
timepoint,  clustered  by  function  and  sorted  by  FDR  value  from  most  to  least  significant. 

Table  S2:  KEGG  pathways  identified  for  differentially  expressed  genes  (FDR  <  0.05).  Pathways 
are  listed  for  each  timepoint  and  sorted  by  FDR  value  from  most  to  least  significant. 


Table  S3.  Significant  GO  clusters  identified  from  the  168h  time  point  that  were  obtained  from 
differentially  expressed  genes  that  were  determined  to  contribute  the  largest  variance  at  that  time 
point. 


Text  SI:  Transcriptional  Signatures  of  Inflammation  and  Immune  System  Detected  in  Early 
Period 

Pro-inflammatory  cytokines  and  chemotactic  protein  members:  IL-ip,  IL-6,  IL-18,  Tnfa,  Osm, 
Ccl2,  Ccr2,  Cxcr2,  Trem2,  S100a8,  FDR=2.21e-22  -  G0:0034097,  response  to  cytokine 

Anti-inflammatory  genes:  Socs3,  IL-lm,  IL-4ra,  IL-lOra,  IL-13ral,  FDR=4.31e-10  - 
G0:0050728,  negative  regulation  of  inflammatory  response 

Invading  immune  cell  genes:  Cd68,  Ly6c,  Cdl4,  Cdl63,  Mrcl,  Ptprc,  Cd24a,  Cd63,  Itgam, 
Clec7a,  Ltf,  Irf7,  Hp,  FDR=7.6e-7  -  G0:0050900,  leukocyte  migration  involved  in  inflammatory 
response 

Cytokines  and  chemokines:  Yml,  S100a9,  Sppl,  Ccrl,  Ccl3,  Cell  1 ,  Lcn2,  Cxcl5,  Ccll2, 
FDR=3.8e-44  -  G0:0001816,  cytokine  production,  FDR=9.05e-14  -  G0:0032602,  chemokine 
production 

Signaling  lymphocytic  activation  molecule  family:  Slamf7,  Slamf8,  Ly9,  FDR=3.44e-16  - 
G0:0046649,  lymphocyte  activation 

Pro-apoptotic  loci:  Caspl,  Casp4,  Casp8,  Apafl,  Fas,  FDR=1.40e-9  -  G0:2001235,  positive 
regulation  of  apoptotic  signaling  pathway 

Anti-apoptotic  loci:  Xiap,  Birc3,  Birc5,  Api5,  Bax,  Mcll,  Bcl211,  FDR=2.17e-23  -  G0:0043066, 
negative  regulation  of  apoptotic  process 

Text  S2:  Transcriptional  Signatures  Associated  With  Injured  Muscle  Tissue  Extracellular  Matrix 

Structural  components  of  extra-cellular  matrix  (ECM)  remodeling:  Collal,  Colla2,  Col3al, 
Col5al,  Col5a2,  Col5a3,  Col6al,  Col6a2,  Col6a3,  Eln,  FDR=4.19e-lll  -  G0:0016043,  cellular 
component  organization 

Connective  tissue:  Fbnl,  Ctgf,  Hasl,  Hbegf,  Fnl,  Sparc,  Tnc,  FDR=6.65e-14  -  G0:0061448, 
connective  tissue  development 

ECM  cross-linking  enzymes:  Lox,  Loxll,  Loxl2,  Plod3,  Leprel2,  P4hb 

Small  leucine-rich  proteoglycans:  Den,  Bgn,  Kera,  Vcan,  Lum,  Ogn,  FDR=3.69e-28  - 
GO: 003 62 1 1,  protein  modification  process 

Adhesive  glycoproteins:  Pcolce,  Itbpbpl,  Lgalsl,  Dpt,  Postn,  Emilin2,  Efemp2,  Igfbp7, 
FDR=3.36e-16  -  G0:0043062,  extracellular  structure  organization 

Tgf-(3  superfamily:  Tgf-pi,  Tgf-Prl,  Tgf-Pr2,  Tgf-pi 

Fonnation  and  stabilization  of  blood  vessels  and  arterial  repair:  Ang,  Angpt2,  Angptll,  Angptl4, 
Angptl7,  Cthrcl,  Aifl,  FDR=3.12e-10  -  G0:0001525,  angiogenesis 

Contractile  fibers:  ActB,  Actg,  Capg,  Arpc3,  Cotll,  Tlnl,  Tagln,  Vcl,  RhoA,  FDR=7.8e-5  - 
G0:0030833,  regulate  actin  filament  polymerization,  FDR=1.5e-3  -  G0:0032273,  positive 
regulation  of  protein  polymerization 

Integrin-associated  genes:  Itga5,  Itga7,  Itgam,  Itgav,  Itgbl,  Itgb2,  FDR=3.34e-5  -  G0:0007229, 
integrin-mediated  signaling  pathway 

Mmps,  Adams  and  Timps:  Mmp2,  Mmp3,  Mmp8,  Mmpl4,  Mmpl9,  Adam8,  Adaml5,  Adaml9, 
Timpl,  Timp2,  FDR=6.83e-8  -  G0:0048771,  tissue  remodeling 


Text  S3:  Injured  Muscle  Tissue  Changes  from  Pro-  to  Anti-Inflammatory  States  Several  Days 
After  Injury 

Innate  immunity  and  microbial  recognition:  Tlrl,  Tlr7,  Tlr8,  FDR=0.003  -  G0:0034121, 
regulation  of  toll-like  receptor  signaling  pathway 

Antigen  presentation  and  hematopoietic  activation  and  proliferation:  Ptprc,  IL-21r,  IL-17ra, 
Vavl,  Lyz2,  Csflr,  CsfZrb,  CsfZra,  I13ra,  I15ra,  Ifi204,  Aifl,  FDR=1.4e-8  -  G0:0042110,  T  cell 
activation,  FDR=8.18e-19  -  G0:0030097,  Haematopoiesis 

Anti-inflammatory  cytokines:  IL-lrll,  Il-4ra,  IL-lOra,  Msrl,  IL-13ral,  S100a9,  FDR=2.9e-10  - 
KEGG  mmu04630,  cytokine-cytokine  receptor  interaction,  FDR=9e-3  -  G0:0010934, 
macrophage  cytokine  production 

Phagocytic  and  complement  cascade  genes:  Cyba,  Ncfl,  Ncf2,  Ncf4,  Clqa,  Clqb,  Clqc,  Clra, 
Cls,  FDR=3.37e-8  -  G0:0006956,  complement  activation,  FDR=7e-6  -  G0:0006909, 
phagocytosis 

Notch  signaling:  Dill,  Notch2,  Myc,  Rbpj,  Cdknla,  FDR  =  3.69e-18,  G0:0002684,  positive 
regulation  of  immune  system  process 

Bone  morphogenetic  proteins:  Bmpl,  Bmp4,  FDR  =  0.08,  G0:0005 125,  cytokine  activity 
Actively  proliferating  cells:  Cdkl,  FDR  =  9.63e-9,  G0:0007049,  cell  cycle 
Text  S4:  Activation  of  Muscle  Repair  Machinery 

Satellite  cell  markers:  CD34,  Itga7,  Itgbl,  Sdc3,  Sdc4,  Sdcbp,  Cavl,  Cxcr4,  Cdhl5,  Ly6a, 
FDR=4.933e-5  -  G0:0007519,  muscle  tissue  development 

Muscle-fusion  genes:  Tmem8c,  Cdhl5,  Capn2,  Capnsl,  Myof,  Cav3,  Itgbl,  Vcaml,  Itga4, 
FDR=2.5e-6  -  G0:0098602,  single  organism  cell  adhesion 

Basement  membrane  genes:  Col4al,  Col4a2,  Dysf,  Trim72,  Lamc2,  Cav3,  Ahnak,  Anxal, 
Anxa2,  Anxa5,  FDR=4.78e-9  -  G0:0061024,  membrane  organization 

Muscle-specific  actins,  myosins,  troponins,  and  tropomyosins:  Actl9,  Clml,  Myh3,  Myl4, 
Myl6b,  Caldl,  Tpm4,  Tnnil,  Tnni3,  Tnntl,  Tnnt2,  Tnncl,  Tpm3,  Des,  Vim 

Transcription  factors:  Itgb2,  Rrad,  Fos,  JunB,  FosB,  Pbx2,  Pbx3,  Dusp5,  MyC,  BMyc,  Mt3, 
Egrl,  Egr2,  Sgms2,  Runxl,  Grn,  Pdlim3,  Csrp3,  Usfl,  Arid5b,  Naca,  Ankrdl,  Ankrd2,  Atf3, 
Atf4,  Atf6,  Xbpl,  Tead4,  Nfe212,  FDR=2.46e-10  -  G0:0006351,  transcription,  DNA-templated 

Text  S5:  Late  Transcriptional  Programs  Activated  in  Response  to  Traumatic  LLMI 

Chemotaxis  genes:  Ccl6,  Ccl8,  Ccl9,  Cxcll6,  Fcerlg,  Fcgr3,  FDR=1.9e-5 

Mast  cell  activation:  Fyb,  Fcerlg,  Fcgr3,  Fcgr2b,  Lat2,  Lcp2,  FDR=3.33e-5  -  G0:0045576,  mast 
cell  activation 

Fc  receptors:  Fcerlg,  Fcgrl,  Fcgr3,  Fcgr2b,  Clec7a,  Colecl2,  Hck,  Sirpa,  Slcllal,  Vavl, 
FDR=0.035  for  KEGG  mmu04664,  Fc  epsilon  RI  signaling  pathway. 

Lectins:  Glee  12a,  Clec4al,  Clec4a2,  Clec4a3,  Clec4n,  Clec7a,  CleclOa. 

Collagen  fibril  organization:  Adamtsl4,  Col3al,  Col5al,  Col5a2,  Dpt,  Lox,  and  Anxa2,  FDR  = 
le-7) 


ECM  receptors:  VLA  proteins  (a2,  a5,  a8,  a9,  alO,  all),  cytoadhesins  (P3,  All-P),  and 
proteoglycans  (Cd44,  Sdc3,  Sv2),  KEGG  ECM-receptor  interaction  pathway,  FDR  =  5.7e-7 

Angiogenesis  genes:  Col4al,  Col4a2,  Robo4,  Vashl,  Sema5a,  Soxl8,  G0:0001525,  FDR  = 
3.24e-14 

Muscle  contraction  and  muscle  system  processes:  Cacnals,  Myoml,  Myh2,  Myh4,  Myh7,  Nosl, 
Kcnmal,  Ryrl,  Trim63,  FDR=6.7e-4  and  1.6e-3,  respectively. 

Ossification  regulation:  Gpnmb,  Bglap,  Sparc,  Sppl,  Runx2,  Bmpl,  Bmp5,  Smadl 


a) 


Fig.  SI.  Immuno-histological  analysis  of  muscle  after  cryo-injury.  a)  Quantitative  analysis  of 
immuno-histological  images  revealed  a  gradual  increase  in  tissue  damage  until  48  hours, 
indicating  secondary  damage  to  the  muscle  tissue  occurred  after  the  initial  cryo-injury.  b) 
Representative  histological  images  of  the  injured  TA  stained  with  Evans  Blue  Dye  (EBD  -  red) 
and  DNA  (blue)  after  several  time  points,  c)  Representative  immuno-histological  images  at  48 
hours  and  672  hours  after  injury,  whereby  the  green  is  stained  for  actin  and  red  is  stained  for 
evans  blue  dye. 
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Fig.  S2.  a)  Example  of  reproducibility  from  RNA-Seq  datasets.  Biological  replicates  of  in  vivo 
isolated  RNA  samples  from  24  hours  after  injury.  Sequencing  libraries  were  prepared  from  the 
isolated  RNA  and  a  strong  correlation  (R2  =  0.95)  was  observed  for  two  representative  examples 
b)  Validation  of  RNA-Seq  data  with  qPCR  comparing  the  injured  versus  uninjured  samples  for 
the  3  hour  time  point,  c)  Validation  of  RNA-Seq  data  with  qPCR  comparing  the  injured  versus 
uninjured  samples  for  the  10  hour  time  point,  d)  Validation  of  RNA-Seq  data  with  qPCR 
comparing  the  injured  versus  uninjured  samples  for  the  168  hr  time  point. 
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Fig.  S3.  Number  of  differentially  expressed  genes 
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b) 


Fig.  S4.  Gene  Sets  Activated  Early  After  Injury  (3-24  hrs)  a)  Enriched  KEGG  pathways  from 
differentially  expressed  genes  for  the  early  time  points  (3-24  hrs).  The  size  of  the  circle 
corresponds  to  the  number  of  significant  genes  with  each  enriched  pathway,  b)  Heatmap  of  genes 
associated  metabolism  expressed  in  the  early  time  period,  c)  Gene  expression  profile  of  Arginase 
1  (Argl),  which  is  a  molecule  emitted  by  type  2  macrophages  to  balance  inflammatory  reactive- 
oxygen  species  present  in  the  injured  tissue.  Red-Injured,  Blue-Uninjured. 
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Fig.  S5.  Examples  of  alternative  splicing  in  immune  network  of  genes  as  determined  by 
Cufflinks. 
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Fig.  S6.  Genes  associated  with  hematopoietic  activation,  proliferation  and  anti-inflammatory 
cytokines  activated  in  early  &  middle  periods. 
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Fig.  S7.  Gene  sets  associated  with  extra-cellular  matrix  remodeling  &  satellite  cell  activation 
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Fig.  S8.  a)  Enriched  KEGG  pathways  from  differentially  expressed  genes  for  the  late  time 
points.  The  size  of  the  circle  corresponds  to  the  number  of  significant  genes  with  each  enriched 
pathway,  b)  Heatmaps  of  gene  sets  activated  in  late  period  categorized  by  function. 
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Fig.  S9.  Area  plots  of  FDR  values  for  Gene  Annotation  (GO)  clusters  over  time.  GO  terms  were 
separated  into  major  functional  categories,  and  the  IQR  was  computed  for  the  terms  in  each 
group.  The  grey  lines  indicate  the  25th  and  75th  percentiles,  red  lines  denote  the  25%  of  gene 
FDR  values,  the  blue  line  denotes  the  median  FDR  value,  and  the  green  line  denotes  the  75%  for 


FDR. 


Timepoint  Clustering  by  PCA  Components  1-2 


Fig.  S10.  Pathway- level  scoring  results  for  test  datasets  from  PCA  and  the  timepoint  signatures 
methods,  a)  Principal  component  analysis,  b)  Normalized  time  point  signatures  method  applied 
to  mean  pathway  expression  values.  The  top  panel  represents  the  timepoint  signatures  obtained 
for  the  uninjured  control  datasets  for  6  test  samples  and  a  truth  training  sample.  The  green  curve 
indicates  the  control  sample  that  was  misclassified  by  the  algorithm  at  the  pathway  level.  The 
next  panels  are  derived  from  time  point  signatures  for  the  9  time  points  (blue  lines)  as  well  as  6 
test  samples  (red  and  green  lines). 


Fig.  Sll.  Differentially  expressed  gene  pathways  in  injured  samples  over  time. 


Table  S2.A  KEGG  pathways  significant  at  3  hours 

Category  Term 


KEGG_PATHWAY  mmu04060:Cytokine-cytokine  receptor  interaction 


KEGG_PATHWAY  mmu04512:ECM-receptor  interaction 

KEGG_PATHWAY  mmu04610:Complement  and  coagulation  cascades 


KEGG_PATHWAY  mmu04510:Focal  adhesion 


KEGG_PATHWAY  mmu04210:Apoptosis 


KEGG_PATHWAY  mmu04062:Chemokine  signaling  pathway 


KEGG_PATHWAY  mmu04010:MAPK  signaling  pathway 


Count  % 


49 


25 


22 


37 


22 


PValue 


4.52  2.96E-10 


2.31  4.36E-09 


2.03  7.69E-08 


3.41  4.52E-07 


2.03  1.18E-06 


OSMR,  CXCR2,  IL17RA,  TGFB2,  CXCR4,  IL1RAP,  CSF3R,  CSF2RB, 
IL1B,  FAS,  IL13RA1,  CSF2RA,  EGFR,  LTBR,  INHBB,  OSM, 
TNFRSF10B,  PPBP,  CCR2,  PDGFRA,  CXCL1,  CSF3,  IL1R2,  CCL3, 
IL1R1,  CCL2,  CSF2RB2,  TNFRSF12A,  CSF1,  CCR1,  IL4RA,  CXCL2, 
CXCL9,  CCL9,  PF4,  CNTFR,  CCL7,  CCL6,  CCL24,  TNFRSF1A, 
TNFRSF1B,  IL10RA,  BMP2,  TGFBR1,  TGFBR2,  CCL11,  CXCL14, 
CXCL16,  BMPR1B 

TNC,  COL3A1,  SDC4,  SDC2,  CD44,  COMP,  COL6A2,  COL6A1, 
COL11A2,  THBS1,  THBS2,  THBS3,  THBS4,  SPP1,  FN1,  TNXB, 
COL4A1,  ITGA3,  ITGA4,  COL5A2,  LAMA4,  ITGA5,  COL1A2, 
LAMC2,  COL1A1 

PLAT,  C7,  C5AR1,  MASP1,  C3,  C4B,  CFB,  F13A1,  SERPING1,  CIS, 
C1QC,  PLAUR,  C1QA,  C1QB,  C1RA,  THBD,  F3,  SERPINE1,  CFH, 
CFD,  PROS1,  PLAU 

MYL2,  PGF,  TNC,  COL3A1,  RAC2,  COMP,  COL6A2,  COL6A1, 
PIK3R5,  ZYX,  PAK1,  COL11A2,  THBS1,  THBS2,  THBS3,  FN1, 

SPP1,  THBS4,  ACTB,  EGFR,  COL4A1,  TNXB,  ACTN1,  ITGA3, 
ITGA4,  FLNC,  BIRC3,  FLNB,  COL5A2,  LAMA4,  ITGA5,  FYN, 
PDGFRA,  COL1A2,  RAP1B,  LAMC2,  COL1A1 
CFLAR,  IL1R1,  CSF2RB2,  RELA,  NFKBIA,  NFKB1,  BIRC3, 

PRKAR2B,  TNFRSF1A,  IRAK3,  TNFRSF10B,  MYD88,  RIPK1, 
IL1RAP,  CASP8,  CASP12,  IL1B,  CSF2RB,  PIK3R5,  FAS,  APAF1, 
MAP3K14 


List  Total  Pop  Hits  Pop  Total  Fold  Enrichment  Bonferroni  Benjamini  FDR 


434  244  5738  2.66  4.88E-08  4.88E-08  3.59E-07 


434  83  5738  3.98  7.20E-07  3.60E-07  5.30E-06 


434  75  5738  3.88  1.27E-05  4.23E-06  9.33E-05 


434  198  5738  2.47  7.45E-05  1.86E-05  5.48E-04 


434  87  5738  3.34  1.95E-04  3.89E-05  1.43E-03 


CXCL1,  ADCY4,  CCL3,  CCL2,  GNAI3,  CCR1,  CXCL2,  CCL9,  CXCL9, 

NFKBIA,  NFKB1,  CXCR2,  PF4,  GNG11,  CCL7,  CCL6,  CCL24,  RAC2, 

TIAM1,  PTK2B,  CXCR4,  PIK3R5,  PAK1,  LYN,  RELA,  STAT3, 

33  3.04  4.28E-06  STAT2,  CCL11,  PPBP,  CXCL14,  CXCL16,  CCR2,  RAP  IB  434  182  5738  2.40  7.05E-04  1.18E-04  5.19E-03 

IL1R2,  IL1R1,  FGFR3,  MRAS,  DUSP10,  HSPA1A,  NFKB1, 

HSPA1B,  TGFB2,  FOS,  TNFRSF1A,  BDNF,  RAC2,  HSPA2, 

MAP3K2,  MAP3K1,  JUND,  MAP3K8,  IL1B,  FAS,  PAK1,  RAPGEF2, 

MYC,  EGFR,  RELA,  TGFBR1,  TGFBR2,  FLNC,  FLNB,  DUSP5, 

MAP4K4,  DUSP2,  GADD45G,  NTRK2,  PDGFRA,  HSPB1,  RAP1B, 

41  3.78  1.37E-05  GADD45B,  MAP3K14,  GADD45A,  CD14,  PLA2G4E  434  265  5738  2.05  2.26E-03  3.23E-04  1.66E-02 


Table  S2.B  KEGG  pathways  significant  at  10  hours 

Category  Term  Count  % 


KEGG_PATHWAY  mmu00190:0xidative  phosphorylation 


KEGG_PATHWAY  mmu05012:Parkinson's  disease 


KEGG_PATHWAY  mmu05016:Huntington's  disease 


95 


KEGG_PATHWAY  mmu05010:Alzheimer's  disease 


93 


KEGG_PATHWAY  mmu04510:Focal  adhesion  82 

KEGG_PATHWAY  mmu00020:Citrate  cycle  (TCA cycle)  23 


KEGG_PATHWAY  mmu04260:Cardiacmusclecontraction 


37 


KEGG_PATHWAY  mmu04512:ECM-receptor  interaction  38 

KEGG_PATHWAY  mmu00280:Valine,leucineandisoleucinedegradation  25 


KEGG_PATHWAY  mmu04010:MAPKsignalingpathway 


90 


KEGG_PATHWAY  mmu04670:Leukocyte  transendothelial  migration  47 


103 


PValue  Genes 

ATP5D,  UQCRC2,  ATP5E,  COX11,  UQCRC1,  COXIO,  ATP5B,  CYC1,  NDUFAB1,  UQCRFS1,  COX5A, 

COX5B,  UQCRQ,  NDUFS7,  NDUFS6,  NDUFS5,  UQCR10,UQCR11,  NDUFS4,  NDUFS8,  ATP5L,  ATP50, 
NDUFS3,  ATP5H,  NDUFS1,  ATP5K,  ATP5J,  COX15,  NDUFB10,  NDUFC2,  COX4I2,  COX4I1,  ATP6V1H, 
NDUFC1,  NDUFA10,  COX6C,  NDUFA11,  ATP6V1A,  UQCRH,  ATP5C1,  UQCRB,  NDUFB3,  NDUFB4, 
NDUFB5,  NDUFB6,  NDUFB7,  NDUFB8,  NDUFB9,  COX7B,  COX7C,ATP6VlB2,  ATP5G1,  COX7A2L, 
ATP5G3,  NDUFB2,  COX6B1,  NDUFA4,  NDUFA5,  NDUFA2,  ATP5J2,  COX7A2,  NDUFA3,  NDUFA8, 
COX7A1,  COX8B,  NDUFA9,  NDUFA6,  NDUFA7,  ATP5F1,  NDUFA1,  PPA2,SDHA,  NDUFV3, SDHB, 

2.10  8.69E-22  ATP6V0E2,NDUFV1,SDHC,NDUFV2,SDHD,ATP5A1 

ATP5D,  UQCRC2,  ATP5E,SNCAIP,  UQCRC1,  ATP5B,  UBE2G1,  CYC1,  SNCA,  NDUFAB1,  UQCRFS1, 
COX5A,  COX5B,  UQCRQ,  NDUFS7,  NDUFS6,  NDUFS5,  CASP3,  UQCR10,  UQCR11,  NDUFS4,  NDUFS8, 
ATP50,  NDUFS3,  ATP5H,  NDUFS1,  ATP5J,  NDUFB10,  SLC25A4,  SLC25A5,  CYCS,  NDUFC2,  COX4I2, 
COX4I1,  NDUFC1,  NDUFA10,  COX6C,  UQCRH,  UBC,ATP5C1,  UBB,  UQCRB,  NDUFB3,  NDUFB4, 

NDUFB5,  NDUFB6,  NDUFB7,  NDUFB8,  NDUFB9,  COX7B,  C0X7C,ATP5G1,  COX7A2L,  ATP5G3,  NDUFB2, 
COX6B1,  NDUFA4,  NDUFA5,  NDUFA2,  COX7A2,  NDUFA3,  NDUFA8,  COX7A1,  COX8B,  NDUFA9, 
NDUFA6,  NDUFA7,ATP5F1,  UBE2L6,VDAC2,VDAC3,  NDUFA1,  SDHA,  NDUFV3,  SDHB,  NDUFV1,  SDHC, 
2.10  6.54E-21  NDUFV2,  SDHD,  APAF1,  ATP5A1 

ATP5D,  UQCRC2,  ATP5E,  CLTB,  UQCRC1,  ATP5B,  CYC1,  NDUFAB1,  REST,  CLTC,  COX5A,  UQCRFS1, 
COX5B,  UQCRQ,  NDUFS7,  GPX1,  NDUFS6,  AP2B1,  NDUFS5,  CASP3,  UQCR10,  UQCRH,  NDUFS4, 

PLCB4,  CASP8,  NDUFS8,  CREB3L2,  ATP50,  NDUFS3,  ATP5H,  NDUFS1,  ATP5J,  NDUFB10,  SLC25A4, 
SLC25A5, TAF4B,  CYCS,  NDUFC2,  COX4I2,  COX4I1,  NDUFC1,  NDUFA10,  COX6C,  UQCRH,  ATP5C1, 
UQCRB,  NDUFB3,  NDUFB4,  NDUFB5,  NDUFB6,  NDUFB7,  NDUFB8,  NDUFB9,  COX7B,  C0X7C,ATP5G1, 
COX7A2L,  ATP5G3,  NDU FB2,  POLR2A,TFAM,COX6Bl,TGM2,TRP53,NDUFA4,NDUFA5,NDUFA2, 
COX7A2,  NDUFA3,  NDUFA8,  COX7A1,  COX8B,  HTT,  NDUFA9,  DNALC1,  NDUFA6,  CREBBP,  NDUFA7, 
ATP5F1,  VDAC2,  VDAC3,  NDUFA1,  SOD2,  NDUFV3,  SDHA,  SDHB,  SP1,  BBC3,  NDUFV1,  SDHC,  BAX, 

2.49  1.19E-18  NDUFV2,  SDHD,APAF1,  ATP5A1 


ListTotal  Pop  Hits  PopTotal  Fold  Enrichment  Bonferroni  Benjamini  FDR 


1276  133 


5738 


2.70  1.24E-18  6.18E-19  8.12E-18 


1276  183 


5738 


2.33  2.24E-16  7.48E-17  1.47E-15 


ATP5D,  UQCRC2,  ATP5E,  UQCRC1,  ATP5B,  SNCA,  CYC1,  NDUFAB1,  UQCRFS1,  COX5A,  COX5B,  UQCRQ, 
NDUFS7,  NDUFS6,  NDUFS5,  CASP3,  UQCR10,  UQCRH,  NDUFS4,  PLCB4,  CASP8,  NDUFS8,  ATP50,  FAS, 
NDUFS3,  ATP5H,  NDUFS1,  ATP5J,  NDUFB10,  ADAM10,  CYCS,  NDUFC2,  COX4I2,  COX4I1,  NDUFC1, 
NDUFA10,  CDK5,  COX6C,  UQCRH,  BACE1,  ATP5C1,  EIF2AK3,  UQCRB,  BID,  NDUFB3,  NDUFB4,  NDUFB5, 
CDK5R1,  NDUFB6,  NDUFB7,  NDUFB8,  NDUFB9,  COX7B,COX7C,ATP5Gl,COX7A2L,ATP5G3,NDUFB2, 
TNFRSF1A,  COX6B1,  GAPDH,  NDUFA4,  LPL,  NDUFA5,  NDUFA2,  COX7A2,  NDUFA3,  NOS1,  NDUFA8, 
COX7A1,  COX8B,  NDUFA9,  NDUFA6,  NDUFA7,  ATP5F1,  CAPN2,  NDUFA1,  ITPR2,  ATF6,  SDHA,  NDUFV3, 
2.44  1.09E-17  SDHB,  LRP1,  ATP2A2,  ATP2A3,  NDUFV1,  SDHC,  NDUFV2,  CASP12,  SDHD,  APAF1,  ATP5A1,  CALM 2 

PDGFB,  PDGFA,  PGF,  BCAR1,  CHAD,  VCL,  ACTG1,  ARHGAP5,  PAK2,  PAK4,  SHC1,  ZYX,  COL11A2, 
COL11A1,  PRKCA,  EGFR,  PARVG,  ACTN4,  BRAF,  ROCK1,  ROCK2,  ACTN1,  ACTN2,  FLNC,  FLNB,  FLNA, 
VASP,  VEGFB,  CCND1,  CCND2,  JUN,  VEGFA,  COL1A2,  PDGFRA,  LAMC2,  LAMC1,  COL1A1,  CAV3,  IBSP, 
MYL2,  TNC,  COL3A1,  COL2A1,  ITGB1,  PXN,  MYL9,  IGF1R,  DOCK1,  LAMB2,  RAC2,  COMP,  ITGB7,  ITGB6, 
COL6A2,  PPP1R12A,  COL6A1,  PIK3R5,THBS1,THBS2,  PIK3R1,THBS3,  FN1,THBS4,  SPP1,  PIK3R2, 

ACTB,  COL4A2,  COL4A1,  FLT1,  TNXB,  MET,  BIRC3,  CAPN2,  COL5A2,  VAV1,  KDR,  VWF,  LAMA4,  ITGA5, 
2.15  1.46E-09  GRLF1,  RAP1B,  DIAP1 

DLST,  AC02,  SUCLG1,  CS,  IDH3B,  ACLY,  DLAT,  OGDH,  PDHB,  IDH3A,  SDHA,  SDHB,  IDH3G,  SDHC,  DLD, 
0.60  8.98E-09  SDHD,  IDH2,  IDH1,  PDHA1,  FH1,  SUCLA2,  MDH2,  MDH1 

UQCRC2,  ATP1B1,  UQCRC1,  MYL2,  ATP1B3,TNNC1,  MYL3,  ATP1B4,  CYC1,C0X7B,  COX7C,  CACNB1, 
COX7A2L,  UQCRFS1,  COX5A,  COX5B,  UQCRQ,  TPM4.TPM3,  UQCR10,  UQCRll,  COX6B1,  ACTC1, 
SLC8A1,  COX7A2,  COX8B,  COX7A1,  COX4I2,  COX4I1,  MYH7,  ATP1A1,  COX6C,  CACNA2D4,  ATP2A2, 
0.97  2.11E-06  UQCRH,  UQCRB,  SLC9A1 


1276  182 


5738 


1276  198  5738 

1276  31  5738 

1276  78  5738 


2.30  2.07E-15  5.17E-16  1.36E-14 


1.86  2.76E-07  5.53E-08  1.82E-06 

3.34  1.70E-06  2.83E-07  1.11E-05 

2.13  3.99E-04  5.70E-05  0.002622723 


IBSP, TNC,  NPNT,  COL3A1,  COL2A1,  SDC4,  ITGB1,  SDC2,  CHAD,  LAMB2,  CD44,  COMP,  ITGB7,  ITGB6, 
COL6A2,  COL6A1,  AGRN,  COL11A2,  THBS1,  COL11A1,  THBS2,  THBS3,  FN1,  SPP1,  THBS4,  COL4A2, 

1.00  4.27E-06  COL4A1,  TNXB,  COL5A2,  VWF,  LAMA4,  SDC1,  CD36,  ITGA5,  COL1A2,  LAMC2,  LAMC1,  COL1A1 

BCAT2,  EHHADH,  ECHS1,  ACAT1,  HADHA,  HIBADH,  AUH,  HADHB,  MCCC2,  MUT,  MCEE,  OXCT1, 

0.66  7.94E-06  MCCC1,  ACAD8,  HADH,  HMGCL,  BCKDHA,  ACAA2,  ACADM,  ACADS,  DBT,  DLD,  AOX1,  HIBCH,  PCCA 

FGF7,  PDGFB,  PDGFA, TGFB3,  NFKB2.TGFB1,  PRKX,TGFB2,  MAP3K6,  FOS,  CASP3,  MAP3K5,  PAK2, 
B230120H23RIK,  MAP3K8,  FAS,  RAPGEF2,  FGF2,  MYC,  CHUK,  MAP2K6,  PRKCA,  EGFR,  BRAF,  RELA, 
RELB,  PTPRR,  ECSIT,  FLNC,  FLNB,  FLNA,  STK3,  MAP4K3,  MAP4K4,  RASGRF2,  ARRB2,  RRAS2,  JUN, 
GADD45G,  PDGFRA,  HSPB1,  MAPK7,  MAP3K14, GADD45B,  GADD45A,  IL1R2,  IL1R1,  FGFR3,  MKNK2, 
DUSP10,  MAP4K2,  CACNB1,  HSPA1A,  GNG12,  HSPA1B,TNFRSF1A,  KRAS,  HSPA2,  RASGRP3,  RAC2, 
ELK4,  MAP3K2,  PLA2G12A,  JUND,  RASGRP2,  HSPA8,  RASA2,  PTPN7, TRP53, TGFBR1, TGFBR2,  NR4A1, 
MAPK11,TAB1,  CACNA2D4,  RPS6KA5,  DUSP5,  NRAS,  DUSP4,  ATF4,  PLA2G4A,  RPS6KA3,  DUSP2, 

2.36  9.02E-06  RPS6KA4,  DUSP1,  RPS6KA2,  NTRK2,  RAP1B,  PLA2G4B,  CD14,  DUSP7 

ITGAL,  GNAI3,  MYL2,  GNAI2,  MMP9,  BCAR1,  CLDN5,  ITGB2,  ITGB1,  MMP2,  PXN,  ITGAM,  VCL,  MYL9, 
ACTG1,  ARHGAP5,  EZR,  RAC2,  PTK2B,  PIK3R5,  MSN,  PIK3R1,  MLLT4,  PIK3R2,  ACTB,  PRKCA,  ICAM1, 
ACTN4,  ROCK1,  NCF2,  NCF1,  ROCK2,  NCF4,  ACTN1,  ACTN2,  MAPK11,  CTNNA1,  VAV1,  VASP,  THY1, 
1.23  3.16E-05  CYBA,  CYBB,  RASSF5,  PLCG2,  GRLF1,  RAP1B,  JAM2 

PDGFB,  PGF,  PDGFA,  MMP9,  TGFB3,  FOXOl,  MMP2,TGFB1,TGFB2,  SLC2A1,  RALB,  CSF3R,  RARA,  FAS, 
RARB,  EGFR,  RET,  BCR,  BRAF,  RELA,  RXRG, CTNNA1,  VEGFB,  JUP,HIF1A,  JUN,  VEGFA,  PDGFRA, 

LAMC2,  LAMC1,  WNT5A,  FGFR3,  SFPI1,  BCL2L1,  ITGB1,  KRAS,  LAMB2,  RAC2,  FH1,  RUNX1,  AXIN2, 
CSF1R,  FN1,  AXIN1,  TRP53,  BMP4,TCF7,  BMP2,  EPAS1,  TGFBR1,  MET,  CREBBP,  TGFBR2,  SMAD3, 
CBLB,  LAMA4,  CDKN1A,  PLCG2,  PPARD,  FGF7,  NFKB2,  FOS,  CCNE1,  CASP3,  CASP8,  FGF2,  MYC, 
CSF2RA,  CHUK,  PRKCA,  AR,  HSP90AA1,  CYCS,  CDK4,  SMO,  CCND1,  MDM2,  WNT11,  WNT9A,  BID,  PML, 
KIT,TPM3,  IGF1R,  PIK3R5,  PIK3R1,  PIK3R2,  CEBPA,  IL6,  COL4A2,  COL4A1,  MSH2,  FZD1,  BIRC3,  FZD5, 
3.93E-05  STAT3,  RALGDS,  NRAS,  HSP90B1,  RASSF5,  BAX,  RASSF1,  JAK1 


1276 

1276 


83 

46 


5738 

5738 


1276  265 


5738 


1276  119 


5738 


1276  323 


2.06  8.06E-04  1.01E-04  0.005297508 

2.44  0.001499078  1.67E-04  0.009856236 


1.53  0.001703641  1.70E-04  0.011202286 


1.78  0.005957607  5.43E-04  0.039252385 


KEGG_PATHWAY  mmu05200:Pathwaysin cancer 


2.70 


5738 


1.43  0.007403353 


6.19E-04  0.04881099 


Category  Term 


KEGG_PATHWAY  mmu00280:Valine,leucineandisoleucine  degradation 


KEGG_PATHWAY  mmu05010:Alzheimer's  disease 

KEGG_PATHWAY  mmu05012:Parkinson's  disease 

KEGG_PATHWAY  mmu04670:Leukocyte  transendothelial  migration 

KEGG_PATHWAY  mmu00190:0xidative  phosphorylation 
KEGG_PATHWAY  mmu00020:Citratecycle(TCAcycle) 

KEGG_PATHWAY  mmu04666:Fcgamma  R-mediated phagocytosis 

KEGG_PATHWAY  mmu04062:Chemokinesignalingpathway 

KEGG_PATHWAY  mmu05016:Huntington's  disease 
KEGG_PATHWAY  mmu00640:Propanoate  metabolism 


28 


Table  S2.D  KEGG  pathways  significant  at  48  hours 

%  PValue  Genes 

BCAT1,  ACADSB,  BCAT2,  ECHS1,  ACAT1,  HIBADH,  HADHA,  AUH,  HADHB,  MCCC2,  MUT,  IVD,  MCCC1, 
MCEE,  OXCT1,  HADH,  BCKDHA,  ACAA2,  ALDH6A1,  ACADM,  ACADS,  BCKDH B,  DBT,  DLD,  ALDH2,  PCCB, 
1.30  2.79E-12  ALDH9A1,  PCCA 


List  Total  Pop  Hits  PopTotal  Fold  Enrichment  Bonferroni  Benjamini  FDR 


820  46  5738  4.26  4.99E-10  4.99E-10  3.43E-09 


UQCRC2,  UQCRC1,  CYC1,  NDUFAB1,  UQCRFS1,  NDUFS7,  CASP3,  NDUFS4,  UQCR11,  CASP9,  APOE, 

NDUFS8,  CASP8,  IL1B,  NDUFS3,  PLCB2,  NDUFS2,  NDUFS1,  ADAM10,  NDUFB10,  CYCS,  NDUFC1, 

NDUFA10,ATP5C1,  UQCRB,  NDUFB3,  BID,  CDK5R1,TM4SF19,  NDUFB8,  NDUFB9,  NDUFB2,TNFRSF1A, 

PPP3CB,  PPP3CA,  LPL,  NDUFA5,  NOS1,  COX7A1,  COX8B,  NDUFA8,  NDUFA9,  NDUFA6,  CACNA1S, 

55  2.56  4.10E-08  SDHA,  NDUFV3,SDHB,  LRP1,  ATP2A2,  NDUFV1,  SDHC,  NDUFV2,  COX6A2,  C0X6A1,APAF1  820  182  5738  2.11  7.33E-06  3.67E-06  5.04E-05 


UQCRC2,  NDUFB3,  UQCRC1,  NDUFB8,  NDUFB9,  CYC1,  UCHL1,  NDUFAB1,  PINK1,  UQCRFS1,  NDUFB2, 
NDUFS7,  CASP3,  UQCR11,  NDUFS4,  CASP9,  NDUFS8,  NDUFS3,  NDUFS2,  NDUFS1,  NDUFA5,  NDUFB10, 
NDUFA8,  COX8B,  COX7A1,  NDUFA9,  SLC25A5,  NDUFA6,  CYCS,  NDUFC1,  NDUFA10,  VDAC3,  VDAC1, 
2.05  6.92E-08  NDUFV3, SDHA, SDHB,  NDUFV1, SDHC,  NDUFV2,  COX6A2, ATP5C1, COX6A1, APAF1,  UQCRB 

ITGAL,  GNAI3,  MYL2,  GNAI2,  MMP9,  ITGB2,  ITGAM,  VCL,  ACTG1,  VCAM1,  EZR,  RAC2,  PTK2B,  CXCR4, 
RHOA,  PIK3R5,  MSN,  RHOH,  ACTB,  PRKCA,  PIK3CG,  ICAM1,  VAV3,  NCF2,  ACTN4,  NCF1,  NCF4,  ACTN1, 
1.81  5.77E-07  MYL12B,  VAV1,  VASP,THY1,  PRKCB,  CYBA,  CYBB,  RASSF5,  MAPK12,  PLCG2,  RAP1B 


820  133 


5738 


820  119 


5738 


2.31  1.24E-05  4.13E-06  8.52E-05 


2.29  1.03E-04  2.58E-05  7.10E-04 


41 


17 


33 


49 


49 


15 


1.53 


2.28 


UQCRC2,  NDUFB3,  UQCRC1,  NDUFB8,  NDUFB9,CYC1,  NDUFAB1,  UQCRFS1,  NDUFB2,  NDUFS7, 
ATP6V0E,  UQCR11,  NDUFS4,  NDUFS8,  NDUFS3,  NDUFS2,  NDUFS1,  ATP5K,TCIRG1,  NDUFA5, 
NDUFB10,  NDUFA8,  COX8B,  COX7A1,  NDUFA9,  NDUFA6,  ATP6V1H,  NDUFC1,  NDUFA10,  NDUFV3, 
8.65E-07  SDHA, SDHB, ATP6V1A, ATP6V0E2,  NDUFV1, SDHC, NDUFV2, COX6A2, ATP5C1,  COX6A1,  UQCRB 

SUCLG1,  CS,  IDH3B,  ACLY,  DLAT,  OGDH,  IDH3A,  SDHA,  SDHB,  IDH3G,  SDHC,  DLD,  IDH2,  PDHA1,  FH1, 
9.45E-07  SUCLA2,  MDH1 

MARCKSL1,  ARPC4,  ARF6,  ARPC5,  DOCK2,  RAC2,  ARPC3,  GAB2,  ARPC2,  PIK3R5,  INPP5D,  PIK3CG, 
PRKCA,  PTPRC,  VAV3,  LYN,  PPAP2C,  NCF1,  LIMK1,  HCK,  SPHK1,  VAV1,  PRKCD,  WAS,  FCGR1,  VASP, 
2.63E-06  PRKCB,  ARPC1B,  PLA2G4A,  FCGR2B,  CFL1,  PLCG2,  MARCKS 

ADCY7,  STAT5B,  CXCR2,  PRKX,  CXCR4,  RHOA,  CSK,  PLCB2,  PIK3CG,  LYN,  NCF1,  PRKCD,  WAS,  ELMOl, 
PRKCB,  CCR5,  PPBP,  ARRB2,  GNB1,  CX3CR1,  CCR2,  GNAI3,  CCL2,  GNAI2,  FGR,  CCR1,  CCL9,  CCL8,  PF4, 
CCL4,  CCL7,  CCL6,  DOCK2,  RAC2,  PTK2B,  SOS2,  PIK3R5,  VAV3,  HCK,  VAV1,  STAT3,  NRAS,  ADCY9, 
9.62E-06  CXCL14,  CXCL13,  CXCL16,  G  NG 10,  JAK2,  RAP1B 

UQCRC2,  CLTA,  UQCRC1,  CYC1,  NDUFAB1,  UQCRFS1,  NDUFS7,  GPX1,  CASP3,  NDUFS4,  UQCR11, 
CASP9,  NDUFS8,  CASP8,  NDUFS3,  PLCB2,  NDUFS2,  NDUFS1,  NDUFB10,SLC25A5,  CYCS,  NDUFC1, 
NDUFA10,ATP5C1,  UQCRB,  NDUFB3,NDUFB8,NDUFB9,NDUFB2,TRP53,NDUFA5,COX7Al,COX8B, 
NDUFA8,  NDUFA9,  NDU FA6,  VDAC3, SOD2,  VDAC1, SDHA,  NDUFV3,SDHB,  BAX, SDHC,  NDUFV1, 
1.13E-05  NDUFV2,  COX6A2,  COX6A1,  APAF1 

ALDH6A1,  LDHB,  ACADM,  SUCLG1,  ECHS1,  ACAT1,  HADHA,  MUT,  ACSS1,  MCEE,  ALDH2,  SUCLA2, 
2.05E-05  PCCB, ALDH9A1,  PCCA 


820  130  5738 

820  31  5738 

820  98  5738 

820  182  5738 


820  183 

820  30 


5738 

5738 


2.21  1.55E-04  3.10E-05  0.001065029 

3.84  1.69E-04  2.82E-05  0.001162609 

2.36  4.71E-04  6.72E-05  0.0032361 

1.88  0.001720856  2.15E-04  0.011840157 


1.87  0.002020134  2.25E-04  0.013901248 

3.50  0.003656954  3.66E-04  0.025184004 


Table  S2.F  KEGG  pathways  significant  at  168  hours 
Category  Term 

KEGG_PATHWAY  mmu04062:Chemokine  signaling  pathway 


PValue  Genes 

CCL2,  FGR,  NCF1,  HCK,  CCL9,  CCL8,  WAS,  CCL7,  CCL17,  PRKCB,  CCL12,  CCR5,  RAC2,  CCR2,  CX3CR1, 
0.65  4.47E-07  PIK3R5,  XCL1,  XCR1 


List  Total  Pop  Hits  Pop  Total  Fold  Enrichment  Bonferroni  Benjamini  FDR 


129  182  5738 


4.78E-05 


4.78E-05  5.02E-04 


KEGG_PATHWAY  mmu04670:Leukocyte  transendothelial  migration  14 

KEGG_PATHWAY  mmu04650:Natural  killer  cell  mediated  cytotoxicity  13 

KEGG_PATHWAY  mmu05322:Systemic  lupus  erythematosus  12 

KEGG_PATHWAY  mmu04060:Cytokine-cytokine  receptor  interaction  18 


1.93E-06  NCF2,  NCF1,  NCF4,  ITGB2,  PRKCB,  VCAM1,  CYBA,  CYBB,  EZR,  RAC2,  PLCG2,  CLDN2,  PIK3R5,  RHOH 
1.46E-05  BID,  FCGR4,  ITGB2,  NCR1,  PRKCB,  FCGR3,  CD48,  RAC2,  PLCG2,  FCER1G,  PIK3R5,  LCP2,  TYROBP 
1.55E-05  C1QA,  C1QB,  C7,  CD86,  C4B,  FCGR4,  H2-DMB1,  HIST1H4I,  C1QC,  FCGR1,  IL10,  FCGR3 

CCL2,  CSF2RB2,  CCL9,  CCL8,  IL10,  CCL7,  CCL17,  CCL12,  CCR5,  IL10RA,  CCR2,  CX3CR1,  CSF2RB,  IL2RG, 
2.48E-05  IL12B,  XCL1,  XCR1,  CSF2RA 


129  119  5738 
129  122  5738 
129  103  5738 

129  244  5738 


5.23  2.06E-04 

4.74  0.001563801 
5.18  0.001652768 


1.03E-04  0.002168346 
5.22E-04  0.016439609 
4.13E-04  0.017375573 


3.28  0.00265068 


5.31E-04  0.027879105 


Table  S2.G  KEGG  pathways  significant  at  336  hours 
Category  Term 

KEGG_PATHWAY  mmu04512:ECM-receptor  interaction 
KEGG_PATHWAY  mmu05416:Viral  myocarditis 
KEGG_PATHWAY  mmu04062:Chemokinesignalingpathway 
KEGG_PATHWAY  mmu04510:Focal  adhesion 


18 

16 

23 

24 


PValue  Genes 

I  BSP,  TNC,  COL3A1,  ITGB5,  ITGA4,  COL5A3,  COL5A2,  SDC2,  COL5A1,  CD44,  COL1A2,  COL6A2,  COL6A1, 
2.36  7.85E-08  C0L1A1,THBS2,THBS3,SPP1,FN1 

H2-K1,  ITGAL,  MYH3,  H2-D1,  H2-DMB1,  ITGB2,  MYH7,  H2-AB1,  MYH8,  CD86,  H2-EB1,  SGCD,  H2-AA, 
2.10  1.24E-05  MYH13,MYH14,ABL2 

PARD3,  CCL2,  ADCY7,  NCF1,  STAT5B,  CXCL9,  CCL9,  CCL8,  GNG11,  PF4,  CCL6,  CCL12,  CCR5,  ADCY9, 
3.01  1.34E-05  TIAM1,  CXCR4,  GNG10,  CXCL16,  GSK3B,  CCR2,  SOS2,  PAK1,  SHC2 

I  BSP,  FLTl.TNC,  COL3A1,  ITGB5,  ITGA4,  COL5A3,  COL5A2,  COL5A1,  GSK3B,  BCL2,  SOS2,  COL6A2, 

3.15  1.64E-05  GRLF1,  COL1A2,  COL6A1,  COL1A1,  PAK1,  EGF,  SHC2,  THBS2,THBS3,  SPP1,  FN1 


ListTotal  Pop  Hits  PopTotal  Fold  Enrichment  Bonferroni  Benjamini  FDR 


254  83  5738 
254  94  5738 
254  182  5738 
254  198  5738 


4.90  1.08E-05 

3.85  0.001713391 

2.85  0.001841954 
2.74  0.002258263 


1.08E-05  9.24E-05 

8.57E-04  0.014623985 
6.14E-04  0.015722211 
5.65E-04  0.019279342 


Table  S2.H  KEGG  pathways  significant  at  504  hours 
Category  Term 

KEGG_PATHWAY  mmu04512:ECM-receptorinteraction 
KEGG_PATHWAY  mmu05416:Viral  myocarditis 
KEGG_PATHWAY  mmu04062:Chemokinesignaling pathway 
KEGG_PATHWAY  mmu04510:Focal  adhesion 


Count  % 
18 


16 

23 

24 


PValue  Genes 

IBSP,TNC,  COL3A1,  ITGB5,  ITGA4,  COL5A3,  COL5A2,  SDC2,  COL5A1,  CD44,  COL1A2,  COL6A2,  COL6A1, 
2.36  7.85E-08  C0L1A1,THBS2,THBS3,SPP1,FN1 

H2-K1,  ITGAL,  MYH3,  H2-D1,  H2-DMB1,  ITGB2,  MYH7,  H2-AB1,  MYH8,  CD86,  H2-EB1,  SGCD,  H2-AA, 
2.10  1.24E-05  MYH13,  MYH14,  ABL2 

PARD3,  CCL2,  ADCY7,  NCF1,  STAT5B,  CXCL9,  CCL9,  CCL8,  GNG11,  PF4,  CCL6,  CCL12,  CCR5,  ADCY9, 
3.01  1.34E-05  TIAM1,  CXCR4,  GNG10,  CXCL16,GSK3B,  CCR2,  SOS2,  PAK1,  SHC2 

IBSP,  FLT1, TNC,  COL3A1,  ITGB5,  ITGA4,  COL5A3,  COL5A2,  COL5A1,  GSK3B,  BCL2,  SOS2,  COL6A2, 

3.15  1.64E-05  GRLF1,  COL1A2,  COL6A1,  COL1A1,  PAK1,  EGF,  SHC2,  THBS2,THBS3,  SPP1,  FN1 


Table  S2.H  KEGG  pathways  significant  at  504  hours 

List  Total  Pop  Hits  Pop  Total  Fold  Enrichment  Bonferroni  Benjamini  FDR 


254  83  5738 
254  94  5738 
254  182  5738 
254  198  5738 


4.90  1.08E-05 

3.85  0.001713391 

2.85  0.001841954 
2.74  0.002258263 


1.08E-05  9.24E-05 

8.57E-04  0.014623985 
6.14E-04  0.015722211 
5.65E-04  0.019279342 


Table  S2.I  KEGG  pathways  significant  at  672  hours 
Category  Term 

KEGG_PATHWAY  mmu04512:ECM-receptor  interaction 


Count  %  PValue  Genes 

I  BSP,  TNC,  COL3A1,  ITGA4,  COL5A2,  SDC2,  COL5A1,  CD44,  COMP,  COL1A2,  COL6A2,  COL6A1, 
17  2.51  7.58E-08  COL1A1,  THBS2,  COL11A1,  THBS3,  FN1 


List  Total  Pop  Hits  Pop  Total  Fold  Enrichment  Bonferroni  Benjamini  FDR 

223  83  5738  5.27  1.08E-05  1.08E-05  8.96E-05 


I  BSP,  CAV2,  CAV1,  MYL2,  TNC,  COL3A1,  BCL2,  COMP,  COL6A2,  COL6A1,  PAK1,  EGF,  COL11A1,  THBS2, 


KEGG_PATHWAY  mmu04510:Focal  adhesion  26  3.85  1.21E-07  THBS3,  FN1,  VAV3,  IGF1,  ITGA4,  COL5A2,  COL5A1,  PDGFRA,  COL1A2,  GRLF1,  COL1A1,  PARVA  223  198  5738 

KEGG_PATHWAY  mmu04610:Complement  and  coagulation  cascades  15  2.22  7.70E-07  C3AR1,  CFB,  F13A1,  SERPING1,  CIS,  C1QC,  C1QA,  C1QB,  C1RA,  CD55,  F3,  CFH,  C2,  CFD,  PROS1  223  75  5738 


3.38  1.72E-05  8.62E-06  1.44E-04 

5.15  1.09E-04  3.64E-05  9.10E-04 


Table  S3.  Significant  GO  clusters  identified  from  the  168h  time  point  that  were  obtained  from 
differentially  expressed  genes  that  were  determined  to  contribute  the  largest  variance  at  that  time  point. 


Cluster  1:  Enrichment  Score:  11.54 

Annotation 

GO  Term 

Gene  Count 

P-value 

FDR 

GOTERM_BP_FAT 

vasculature  development 

85 

5.20E-14 

9.70E-11 

GOTERM_BP_FAT 

blood  vessel  development 

83 

1.00E-13 

1.90E-10 

GOTERM_BP_FAT 

blood  vessel  morphogenesis 

69 

4.30E-12 

8.00E-09 

GOTERM BP FAT 

angiogenesis 

48 

3.00E-09 

5.60E-06 

Cluster  2:  Enrichment  Score:  11.31 

Annotation 

GO  Term 

Gene  Count 

P-value 

FDR 

GOTERM_CC_FAT 

extracellular  matrix 

100 

2.50E-14 

3.60E-11 

SP_PIR_KEYWORD 

extracellular  matrix 

76 

7.40E-14 

1.10E-10 

GOTERM_CC_FAT 

proteinaceous  extracellular  n 

96 

9.40E-14 

1.40E-10 

GOTERM_CC_FAT 

extracellular  region  part 

179 

1.10E-09 

1.60E-06 

GOTERM CC FAT 

extracellular  matrix  part 

37 

1.50E-08 

2.20E-05 

Cluster  3:  Enrichment  Score:  9.38 

Annotation 

GO  Term 

Gene  Count 

P-value 

FDR 

KEGG_PATHWAY 

lysosome 

53 

2.40E-12 

3.00E-09 

SP  PIR  KEYWORD 

lysosome 

52 

5.60E-10 

8.30E-07 

GOTERM  CC  FAT 

lysosome 

60 

1.00E-09 

1.50E-06 

GOTERM  CC  FAT 

lytic  vacuole 

60 

1.30E-09 

1.90E-06 

GOTERM CC FAT 

vacuole 

64 

6.80E-09 

9.90E-06 

Cluster  4:  Enrichment  Score:  7.85 

Annotation 

GO  Term 

Gene  Count 

P-value 

FDR 

GOTERM  BP  FAT 

vesicle-mediated  transport 

121 

2.70E-10 

5.00E-07 

GOTERM  BP  FAT 

membrane  invagination 

61 

2.10E-09 

3.90E-06 

GOTERM  BP  FAT 

endocytosis 

61 

2.10E-09 

3.90E-06 

GOTERM_BP_FAT 

membrane  organization 

76 

3.10E-08 

5.70E-05 

SP PIR KEYWORD 

endocytosis 

31 

1.60E-05 

2.40E-02 

Cluster  5:  Enrichment  Score:  7.18 

Annotation 

GO  Term 

Gene  Count 

P-value 

FDR 

GOTERM  BP  FAT 

actin  filament-based  process 

61 

1.10E-10 

2.00E-07 

GOTERM  BP  FAT 

actin  cytoskeleton  organizati< 

57 

5.60E-10 

1.00E-06 

GOTERM  BP  FAT 

actin  filament  organization 

23 

6.70E-06 

1.20E-02 

GOTERM_BP_FAT 

cytoskeleton  organization 

76 

4.80E-05 

8.80E-02 

